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-large-radius  concave  surface.  It  was  recognized  that  neither  a geometrical 
optics  nor  a whispering  gallery  mode  representation  is  satisfactory  for 
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integral  for  convex  surfaces,  may  be  the  most  effective.  However,  the 
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appropriate  combination  of  rays  and  modes  were  not  at  hand,  nor  were  compre- 
hensive numerical  comparisons  to  assess  the  utility  and  quality  of  various 
field  representations.  Missing  also  were  tractable  results  for  the  case  when 
the  observation  point  approaches  the  source  or  when  the  radius  of  curvature  of 
a circular  contour  tends  to  infinity  (infinite  plane  limit).  This  work  has 
succeeded  in  establishing  a sound  analytical  basis  for  any  of  the  above- 
mentioned  field  formulations. 
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The  findings  in  this  report  are  not  to  he  construed 
as  an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


The  objective  of  the  research  program  under  Grant  No.  DAHC 


O-i-75-G-OlfiZ,  extending  from  1 June  1975  to  31  May  1978,  has  been  to  gain  an 
under  standing  of  the  electromagnetic  propagation  and  diffraction  properties  of 
concave  surfaceshaving  radii  of  curvature  large  compared  to  the  wavelength, 
with  application  to  ground  wave  propagation,  large  reflector  antennas,  and 
mutually  coupled  conformal  antenna  arrays. 

This  basic  objective  lias  been  accomplished.  When  the  program  was 
initiated,  there  was  only  incomplete  and  non-conclusive  information  available 
on  how  the  radiation,  diffraction  and  guiding  of  fields  is  affected  by  a nearby 
large-radius  concave  surface.  It  was  recognized  that  neither  a geometrical 
optics  nor  a whispering  gallery  mode  representation  is  satisfactory  for  descrip- 
tion of  the  field,  and  that  a mixture  of  rays  and  modes,  or  a more  complicated 
formulation  involving  a canonical  integral  analogous  to  the  Fock  integral  for 
convex  surfaces,  may  be  the  most  effective.  However,  the  foundations  for  these 
conjectures  were  lacking;  accurate  criteria  for  the  appropriate  combination  of 
rays  and  modes  were  not  at  hand,  nor  were  comprehensive  numerical  comparisons 
to  assess  the  utility  and  quality  of  various  field  representations.  Missing  also 
were  tractable  results  for  the  case  when  the  observation  point  approaches  the 
source  or  when  the  radius  of  curvature  of  a circular  contour  tends  to  infinity 
(infinite  plane  limit). 

Our  work  has  succeeded  in  establishing  a sound  analytical  basis  for 
any  of  the  above-mentioned  field  formulations.  Most  important,  we  have  shown 
that  a mixture  of  ray-optical  and  whispering  gallery  mode  fields  alone,  chosen 
in  correct  proportion,  is  indeed  capable  of  providing  accurate  field  values  on 
the  concave  surface  (so  far  taken  to  bo  a circular  cylinder).  Exhaustive  numeri- 
cal calculations  have  provided  reliable  comparisons  of  different  field  representa- 

• • • ^ 

tions  and  have  permitted  a critical  evaluation  of  their  accuracy  and  sensitivity 
to  error  or  to  parameter  changes.  The  results  produced  under  this 


contract  are  the  most  complete  yet  available,  and  they  have  yielded  basic 
insights  into  the  physical  mechanisms  that  establish  the  concave  surface  field. 

Two  manuscripts,  accepted  for  publication  in  IEEE  Transactions  on  Antennas 
and  Propagation,  and  attached  in  the  Appendix,  summarize  these  accomplishments 
which  are  therefore  not  detailed  here: 

1.  "High  Frequency  Fields  Excited  by  a Line  Source  Located  on  a Perfectly 
Conducting  Concave  Cylindrical  Surface",  by  T.  Ishihara,  L.  B.  Felsen  and 

A.  Green,  and 

2.  "High  Frequency  Fields  Excited  by  a Line  Source  Located  on  a Concave 
Cylindrical  Impedance  Surface",  by  T.  Ishihara  and  L.  B.  Felsen. 

A third  manuscript,  almost  completed,  extends  these  studies  to  the  three- 
dimensional  case  of  point  source  excitation. 

The  following  papers  were  given  at  technical  meetings: 

"Excitation  of  Large  Concave  Surfaces",  by  L.B,  Felsen  and  A.  Green, 
IEEE/AP-S  International  Symposium,  University  of  Massachusetts,  Amherst, 

Mass.,  October  10-15,  1976. 

"Ground  Wave  Propagation  in  the  Presence  of  Smooth  Hills  and  Depressions”, 
by  L.B.  Felsen  and  A.  Green,  AGARD  Symposium  on  EM  Propagation  Characteristics 
of  Surface  Materials  and  Interface  Aspects,  held  in  Istanbul,  Turkey,  October  18-19,1976 
and  published  in  the  Symposium  Proceedings. 

"Effects  of  Smooth  Elevations  and  Depressions  on  Ground  Wave  Propagation", 
by  L.B.  Felsen  and  T.  Ishihara,  URSI  Symposium  on  Propagation  in  Non-Ionized 
Media,  held  in  La  Baule,  France,  April  28-May  6,  1977. 

"High-Frequency  Currents  Excited  by  a Line  Source  Located  on  a Concave 
Cylindrical  Surface",  by  L.B.  Felsen  and  T.  Ishihara,  URSI  International  EM  Wave 
Theory  Symposium,  Stanford  University,  Palo  Alto,  California,  June  20-24,  1977. 


"High- Frequency  Fields  Excited  by  a Line  Source  Located  on  a Concave 


Cylindrical  Impedance  Surface",  by  T.  Ishihara  and  L,  B.  Felsen,  1978 
International  URSI/AP-S  Symposium,  Washington,  D.  C.,  May  15-19,  1978. 

"High-Frequency  Behavior  of  Concave  Surfaces"  by  L.  B.  Felsen, 
Invited  paper,  18  URSI  General  Assembly,  Helsinki,  Finland,  August  1978. 
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1.  High  Frequency  Fields  Excited  by  a Line  Source  Located  on  a 
Perfectly  Conducting  Concave  Cylindrical  Surface,  T.  Ishihara, 

L.  B.  Felsen  and  A.  Green. 

2.  High  Frequency  Fields  Excited  by  a Line  Source  Located  on  a 
Concave  Cylindrical  Impedance  Surface,  T.  Ishihara  and  L.  B.  Felsen. 


1 . High  Frequency  Fields  Excited  by  a Line  Source  Located  on  a 
Perfectly  Conducting  Concave  Cylindrical  Surface 


T.  Ishihara,  L.  B.  Felsen  and  A.  Green 
Department  of  Electrical  Engineering 
Polytechnic  Institute  of  New  York 
Farmingdale,  New  York  117  35 


Abstract 


Alternative  representations  are  obtained  for  the  high  frequency  surface 
field  excited  on  a perfectly  conducting  concave  circular  cylinder  by  an 
axial  magnetic  line  current  located  on  the  surface.  Included  are  ray 
optical,  canonical  integral,  whispering  gallery  mode  and  near  field  formu- 
lations, and  various  combinations  of  these.  Asymptotic  evaluations  in 
different  parameter  ranges  lead  to  results  with  varying  accuracy  and 
physical  content.  Their  utility  is  assessed  by  extensive  numerical  cal- 
culations and  comparisons.  Most  intriguing  is  a form  of  the  asymptotic 
solution  that  involves  only  a number  of  geometric  optical  rays  and  a 
number  of  whispering  gallery  modes. 


This  work  was  sponsored  by  the  U.  S.  Army  Research  Office,  Durham, 
North  Carolina,  under  Contract  No.  DAI1C  0-1  -75  -G -0 1 5 2. 
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I.  Introduction 

When  high  frequency  fields  impinge  on  a concave  perfectly  conducting 
surface,  the  induced  surface  currents  can  be  well  approximated  by  their 
physical  optics  values  since  the  entire  surface  is  illuminated.  Phrased 
alternatively,  the  method  of  geometrical  optics  is  adequate  to  describe 
the  surface  field  behavior.  However,  when  actual  sources  or  induced 
equivalent  sources  (for  example,  at  an  edge  termination)  are  located  very 
near  or  on  the  surface,  geometrical  optics  becomes  inapplicable  because 
the  multiply  reflected  ray  fields  (Fig.  1)  have  caustics  that  lie  near  the 
boundary  (Fig.  2).  Therefore,  it  is  necessary  to  account  in  some  other 
manner  for  the  fields  that  would  be  contributed  by  rays  having  undergone 
many  reflections. 

A systematic  study  of  the  problem  can  be  undertaken  on  the  prototype 
configuration  of  a circular  cylinder.  Since  a thorough  understanding  of  the 
axially  independent  two-dimensional  case  is  essential  for  extension  to 
three-dimensional  fields,  attention  is  focused  here  on  the  problem  of 
excitation  by  a magnetic  line  current  located  on  the  concave  surface  and 
oriented  parallel  to  the  cylinder  axis.  Starting  from  a rigorous  formulation 
of  the  Green's  function  problem,  ®ne  may  extract  from  the  aximuthally 
periodic  solution  for  the  closed  surface  a portion*  that  describes  only  the 
propagation  characteristics  between  points  Q and  P in  Fig.  1.  Results 
obtained  therefrom  are  applicable  to  concave  surface  segments  that  are 
open,  and  they  are  in  a form  suitable  for  subsequent  generalization  to 
variable  surface  properties.  This  basic  representation,  given  as  a con- 
tour integral,  is  the  starting  point  for  the  present  study,  which  has  as  its 
goal  the  exploration  of  alternative  field  representations  that  are  useful  for 
calculation  in  various  parameter  ranges.  Included  are  ray-optical,  canonical 
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integral  and  whispering  gallery  mode  (Fig.  3)  solutions,  and  various 
combinations  of  these.  Strong  emphasis  is  placed  on  the  physical  inter- 
pretation of  the  results,  and  on  their  relative  accuracy.  Most  pleasing 
from  a physical  standpoint  is  a field  solution  that  comprises  a mixture  of 
(N+l)  geometric  optical  rays  and  M whispering  gallery  modes,  with  criteria 
provided  for  the  proper  choice  of  N and  M. 

In  the  course  of  the  presentation,  reference  is  made  to  the  sparse 
relevant  literature  on  this  problem'  For  certain  of  the  alternative 

representations,  we  have  performed  extensions  and  refinements  of  previous 
results  but  others  obtained  here,  in  particular,  the  above-mentioned  ray  and 
whispering  gallery  mode  combination,  are  new.  The  conclusions  derived 
from  this  study,  supported  by  extensive  numerical  calculations,  provide  a 
thorough  understanding,  both  physically  and  quantitatively,  of  the  propa- 
gation characteristics  of  perfectly  conducting  concave  cylindrical  surfaces. 
Some  extensions  to  non-circular  and(or)  surface  impedance  boundaries 
have  been  considered  elsewhere  ( as  has  the  case  of  point  source 


excitation  . 
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II.  Alternative  Field  Presentations 

A.  Green's  Function  Formulation 

We  consider  the  prototype  problem  of  line  source  radiation  inside  a 
perfectly  conducting  circular  cylinder  of  radius  a_.  To  make  the  analysis 
relevant  to  the  propagation  along  concave  surfaces  in  general,  it  is  neces- 
sary to  remove  the  azimuthal  periodicity  imposed  on  the  fields  in  the  cylin- 
drical p = (p  , <{>)  geometry,  where  p represents  the  radial  and  <*>  the 
angular  (azimuthal)  coordinate.  This  can  be  done  by  extending  the  range 
of  the  <{>  coordinate  from  its  physical  periodic  domain  0 <_  4>  <_  2 rr  into  an 
infinite  domain  -x<<$>  < x Such  an  extension  implies  that  angularly 
propagating  waves  originating  at  the  source  point  ( p * . 4>' ) are  outgoing 
toward  |<(>|  = x;  this  "angular  radiation  condition"  can  be  realized  by  placing 
at  some  angular  location  away  from  the  source  angle  4>'  a perfect  absorber 
for  angularly  propagating  waves.  Such  an  absorber  has,  however,  the 
uridesired  property  of  generating  diffraction  at  the  radial  coordinate  origin 
p = 0.  Therefore,  when  considering  propagation  phenomena  ascribable 
only  to  the  cylinder  surface,  it  is  desirable  to  remove  the  spurious  dif- 
fraction effects  from  the  total  field  solution. 

The  line  source  Green' s function  in  the  cylindrical  domain  0 < p < a, 

-ec  < <}>  < x , can  be  constructed  by  the  method  of  separation  of  variables 

and  expressed  in  its  most  general  form  as  a contour  integral  involving  the 

two  one -dimensional  characteristic  Green's  functions  g and  g,  for  the 

bP  6<t> 

radial  and  angular  domains,  respectively^.  When  the  Green's  function  is 
represented  in  terms  of  angularly  propagating  waves,  the  cor  ’ms  ponding 
eigenmode  spectrum  in  the  radial  domain  involves  a discrete  (whispering 
gallery  mode)  and  continuous  portion^.  The  latter  accounts  in  part  for 
propagation  phenomena  associated  with  the  portion  of  the  concabc  boundary 
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lying  between  the  source  and  observation  points  (this  is  the  desired  contribu- 
tion). and  in  part  for  the  spuriou".  diffraction  effects  of  the  angular  absorber. 
When  the  spurious  effects  are  removed  and  cognizance  is  taken  of  the  liigh- 
frequency  nature  of  the  analysis,  one  may  show  that  the  relevant  propagation 
characteristics  along  the  portion  of  the  concave  boundary  lying  between 
the  source  point  Q and  the  observation  point  P in  Fig.  1 are  contained 
in  the  partial  Green's  function 


G(p.p'  ) = 2 

i(Trka) 


iv  | <b-  <>' 


H'v(2,(ka)J^(ka) 


where  k is  the  free-space  wavenumber,  the  prime  on  the  cylinder 

functions  denotes  the  derivative  with  respect  to  the  argument,  and  a time 

factor  exp(-iut)  is  implied.  It  has  been  assumed  in  (1)  that  P and  Q 

both  lie  on  the  boundary  (i.  e.  , p = p*  = a)  and  that  the  boundary  condition 

requires  that  the  normal  derivative  of  the  Green's  function  vanishes  at 

p = a;  this  makes  G proportional  to  the  axial  component  of  magnetic  field. 

The  contour  C and  the  singularities  of  the  integrand  in  the  complex 

v -plane  are  shown  in  Figure  4.  Contributions  from  the  pole  singularities 

v arising  from 
m 

J'  (ka)  =0,  m = 1, 2...  . (z) 

. v 

m 

are  found  to  describe  whispering  gallery  modes.  While  (2)  has  an  infinite 
number  of  real  solutions  as  indicated  in  Fig.  4,  only  those  with  Re  0 

represent  spectral  contributions  in  the  angular  transmission  representation 
which  includes  also  a continuous  spectrum. 
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B.  Whispering  Gallery  Mode  and  Continuous  Spectrum  Representation 
The  spectral  representation  comprising  a discrete  (whispering  gallery 
mode)  portion  and  a continuous  portion  may  be  obtained  directly  from  (1)  by 
deforming  C in  Fig.  4 into  a contour  extending  along  the  imaginary  v axis: 


G = T G + R„. 
m = 1 1 

where  is  the  total  number  of  poles  satisfying  Re  vm  > 0 and 


R T 

M1  i(irka)2  -ioo  H'  (2)(ka)J'  (ka) 

V V 


,1  Je-toU-*'|»inh»dVi  ka>>,  (4a) 


is  the  continuous  spectrum.  Following  Wasylkiwskyj,  the  first  integral  may 
be  simplified  to  the  second  form  shown  and  also  to 


-^[^(ks)  - Yo(ks)]  . s=aU-4,'l 


where  H is  the  Struve  function  and  Y the  Neumann  function.  Each 
— o o 

residue  contribution  G^  represents  a whispering  gallery  mode  of  the  form 


Gm  = ka  ' 


i ‘vjo-*’! 


Jvm<k»  Psr  JV<k*>] 


On  use  of  the  uniform  asymptotic  formulas  given  by 


J (ka)  ~ / r — - 

v A/  ka  cos  w 


a4  Ai[  -a]  . v = ka  sin  w 


Jv'(ka)v>nr  a *4  Ai'  [ -o] 


= -|- ka  [cos  w - - w)  sin  w ] , Rew>0 


i 
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which  are  justified  here,  one  finds 


exp  (iks  sin  \v  + in/2) 
, m 

Jm  ku ( tt / 2 - w )cos  \v 

m m 


The  roots  w are  determined  by  solution  of  the  resonance  equation 

, ,,  . . 2 _ 3/2 

cosw  - C tt / 2 - w ) sinw  = — rr — or 
m m m 3ka  m 

where  Ai'(-am)  = 0 , m = l,2,... 


C.  Ray-optical  Representation 

To  generate  a ray-optical  formulation,  one  replaces  H^Z\ka)  and 
J^(ka)  = { II'  ^ \ka)  + \ka)} /Z  by  their  Debye  asymptotic  forms 

Hy  ’2^(ka)~  exp  ika  [ cos  w - ( j-  - w)sin  w ]+_i  ^ J (8) 

Re  w > 0,  ^ka  [ cos  w • ( —•  - w)  sin  w > > 1 (8a 

where  v ? ka  sin  w.  Utilising  the  traveling  wave  expansion 


1 . 2 S „ n _ Hj/V, 

Jv(ka|  IT<2)(ka)  it^O  ’ fr,2,(k.i) 


. I r I < 1 . 


In  conjunction  with  (8),  one  reduces  (1)  to  a series  of  integrals 


where 


G ~ — (*i)n  f exp  I ika  q (w)l  dw 
n n A , r ‘ 'n  J 


(10a) 
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and 

q^(w)  = | £-<*>'  | sin  w + 2(n+l)  { cos  w - ( - w)  sin  w]  . 

The  saddle  points  of  q (sv),  as  obtained  from  dq  /dw  = 0,  are 

n 1' 


it  I <fr  -6 1 1 ~ _ JT 

Wsn  " 2 ' 2(n4  1)  1 s ' 2 


(10b) 


(U) 


and  thus  lie  on  the  real  \v  axis  between  w = 0 and  w = tr/2.  A typical  path 
through  the  nth  saddle  point  is  shown  in  Figure  5.  The  original  path  mapped 
.from  the  v plane  can  be  deformed  into  C\  The  Debye  approximations  in  (8) 
are  valid  in  this  relevant  region  of  the  complex  w-plane.  Use  of  the  conven- 
tional saddle  point  formula  for  evaluation  of  the  integral  in  ( 10a)  yields 


G 

n 


ikD 

n 


(12) 


where 

= 2(n+  1)  a sin  [ | <$>  -<{>  1 1 /2(n+  1)  ] ( 12a) 

This  result  corresponds  to  a direct  or  multiply  reflected  ray  as  depicted 
in  Figure  1.  and  could  have  been  constructed  directly  by  ray-optical 
techniques. 

Although  the  series  in  (10),  with  (12),  formally  contains  ray  contribu- 
tions with  an  arbitrarily  large  number  of  reflections,  these  are  suspect 
since  the  saddle  points,  from  which  they  are  derive^  all  cluster  about 
w = ir/2.  Thus,  the  asymptotic  method,  whereby  each  saddle  point  is  treated 
as  Isolated,  is  inadequate.  Moreover,  (8)  becomes  invalid  as  w — • tt/2, 
thereby  Invalidating  the  simplification  of  the  integrand,  on  which  the  saddle 
point  evaluation  is  based.  It  is  therefore  necessary  to  truncate  the  number 
of  legitimate  ray-optical  terms  at  some  n = N such  that  w^^  sufficiently 
less  than  ir/2. 


I 


D.  Ray  plus  Canonical  Integral  Representation 

In  view  of  the  above  observations,  one  may  employ  instead  of  19 ) the 
partial  expansion 


J1  (ka) 

v 


H'  (2)(ka)  n*=0 


^ (-1)"  r"  + -i 


J' (ka) 
v 


. ,,N+1  N+l 
(-1)  r 


Substitution  of  (13)  into  (1)  yields 


N 


G = V G + R., 
n N 

n=0 


where 


_1 

i(trka)2 


l 


| ‘l'  - 4> ' I 

H'  (2)(ka)J'(ka) 

v v 


N+  1 
r 


(-l)N+1dv 


(13) 


(14) 


(15) 


The  sum  in  (14)  evidently  describes  ray  contributions  having  experienced 
up  to  N reflections.  The  remainder  integral  incorporates  the  cumulative 
effects  of  ray  fields  having  been  reflected  more  than  N times.  It  should  be 
noted  that  when  the  observation  point  approaches  the  source  point  so  that 
| - 4>  * | — • 0 , one  has  vs  — ka  even  for  n = 0,  and  the  ray  representation 
fails  altogether.  When  N = - 1,  no  ray  contributions  are  extracted  and  (1  5) 
agrees  with  ( 1). 

Several  options  are  open  for  dealing  with  R^.  One  possibility  is  to 
explore  under  what  conditions  the  contributing  range  of  the  integrand  is 
localized  so  that  approximation  methods  may  be  used  effectively.  To  this 
end,  the  cylinder  functions  are  first  replaced  by  their  uniform  asymptotic 
approximations  in  terms  of  Airy  functions.  From  the  behavior  of  the  inte- 
grand over  the  contour  C in  the  v-pl^nc  or  the  contour  C*  in  the  w-plane, 

It  may  then  be  shown  that  the  principal  contribution  to  the  integral  arises 
from  the  vicinity  of  v = ka  or  w = n/2  provided  that  the  arc  length  parameter 
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( ka  1 /3S 

r = ( — ) 


( lCa) 


2 • a 

and  the  number  of  rays  are  related  by  the  criterion  (note  that  \vc  designate 


a ray  having  undergone  N reflections  as  the  (N+l)th  ray) 

r Y/T  - 1 < (N+  1)  < r/2 

2(21/3Ar~  ~ 2(ZUi±)' 


(16b) 


where  A = A(ka)  = 0(1)  is  chosen  appropriately.  When  the  left-hand  side  of 
the  inequality  is  negative,  i.  e.  , in  the  region  near  the  source,  no  ray- 
optical  contributions  are  included.  The  first  inequality  in  (16b)  is  required 
to  assure  that  the  integrand  in  (18a)  decays  away  from  t = 0.  The  second 
inequality  in  (16b)  (and  the  inequality  in  (21a))  is  required  to  validate  use  of 
the  Debye  approximations  in  the  first  term  in  (13),  which  yields  the  geometric 
optical  ray  contribution. 

Our  criterion  in  (16b),  the  utility  of  which  is  confirmed  by  the  numerical 
calculations  in  Sec.  Ill,  differs  somewhat  from  that  obtained  by  Babich  and 
Buldyrev^.  The  numerical  factor  A in  the  inequalities  (16b),  and  also  in 
(21a),  is  not  rigidly  fixed.  This  implies  that  the  location  of  the  point  at 
which  a geometric-optical  ray  field  splits  off  from  the  remainder  field 
can  be  varied.  In  each  specific  problem,  the  factor  A should  be  chosen  ( 

such  that  the  discontinuity  in  the  total  field  is  a minimum  when  the  new 
■geometric-optical  wave  emerges.  The  estimated  value  of  A = 0(1)  can  be 
determined  from  a knowledge  of  the  point  of  emergence  of  the  direct  ray. 

This  point  can  be  found  by  comparing  the  numerical  results  of  the  direct 
ray  plus  canonical  integral  representation  here,  or  of  the  direct  ray  plus 
whispering  gallery  mode  representation  (Sec.  E).  with  that  of  the  exact  solution 
(whispering  gallery  mode  and  continuous  spectrum  representation).  In  this 
manner,  one  finds  from  ( 1 6b)  or  (21a)  with  N = 0 that  y = 3.  1 whence 
A = 2 (see  Figs.  (8),  (10)  and  (11)).  This  value  of  A is  then  used  in  (16b) 
to  determine  the  splitting  off  points  of  the  various  reflected  rays.  It  turns 


ir 
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out  that  the  results  are  not  very  sensitive  to  the  choice  of  A because  of 
the  overlap  of  the  curves  for  various  N.  Therefore,  A can  be  determined 
from  the  first  overlap  region  for  N = -l,  N = 0,  . without  the  need  of  per- 
forming the  comparison  with  the  exact  solution. 

Changing  variables  to  t via  v = ka  + (ka/2)1//3t  and  using 


. . 

H'(2)(ka)  ~ - (-r-)2/3  w'  (t).  J’  (ka)  ~ - (-j^)2/3  Ai'l 


where  w^(t)  = Ai(t)  + iBi(tVone  may  write  (15)  as 


N + l eikS  2/3 


2iir  2 ka 


where 


( iN(v) 


eiYt  rw'(t)“|N+l 
£ Ai'-(t)  w^(t)~  [_  w'2(t)  J 


IN(V)  - 


(18a) 


The  contour  in  the  t-plane  is  inferred  from  the  mapping  and  Fig.  4 by 
observing  that  v = ka  corresponds  to  t = 0.  The  integral  in  (18a)  is  ir  a 
canonical  form  analogous  to  the  Fock  integral  for  convex  surfaces.  Tabula- 

Q 

tions  exist  for  selected  ranges  of  y and  N«  When  N = -1,  the  resulting  ex- 
pression in  (18)  is  a valid  representation  for  G in  the  range  of  small  y. 

N = 0 defines  the  domain  of  the  direct  ray,  N = 1 the  domain  of  the  direct 
plus  singly  reflected  ray,  etc. 

E.  Ray  plus  Whispering  Gallery  Mode  Representation 
An  alternative  possibility  is  to  represent  R^  as  a whispering  gallery 
mode  sum  plus  an  explicit  remainder.  If  C is  deformed  into  in  Fig.  4. 

M residues  at  the  poles  v of  the  integrand  lying  between  C and 
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must  be  extracted.  One  finds 


^ = ^ Gm  + RMN 

m=l 


(19) 


where  G is  given  in  (5)  or  (7).  The  remainder  integral  R^j^  is  the 
same  as  except  for  the  replacement  of  C by  C^.  It  may  be  shown 

(see  Appendix)  that  the  remainder  integral  can  be  approximated  by 


RMN 


(20) 


provided  that  observation  poinsts  are  characterized  by  wsn  * ( wj^£  + WM+1^’ 
i.  e.  , for  saddle  points  wg^  well  removed  from  the  pole  singularities 
(see  Section  III  for  a discussion  on  procedure  when  this  condition  is  not 
satisfied).  Thus. 


G 


G 


m 


(21) 


The  number  of  rays,  (N+l),  is  picked  according  to  the  criterion 


N + 1 < 


2(21^3A)1/,Z 


(21a) 


F.  Near  Field  and  Infinite  Plane  Limits 

Except  for  the  whispering  gallery  mode  plus  continuous  spectrum 
representation  in  Sec.  B,  which  is  not  convenient  for  numerical  computation, 
the  other  formulations  cannot  account  for  the  limiting  case  y -»  0 in  (16a), 
which  arises  either  when  the  observation  point  approaches  the  source  point 
(s  -*  0 with  fixed)  or  when  the  radius  of  curvature  becomes  arbitrarily 
large  (a  -*  oo  with  s arbitrary  but  fixed).  The  latter  limit  traces  the  transi- 
tion from  a concave  to  a plane  surface  and,  if  a^  is  allowed  to  become 
negative,  from  concave  to  convex.  We  have  derived  the  transition  functions 
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4 


for  these  cases. 


We  begin  with  G expressed  by  (IS),  with  N = -1,  and  use  the  Wronskian 
relation  for  the  Airy  functions  together  with  Cauchy's  theorem  applied  to 
the  upper  half  of  the  t-plane  to  infer  that,  equivalently. 

--^<¥>2/3  X •**  m * 


Employing  large  argument  expansions,  one  may  show  that 


Aiit]  _ f°  ..  a| 

Ai'(t)  t(l+3j)/2 


+ 0(t'34,/V  - ir  < arg  t < 


where  aQ,  a^,  ...  a^  are  given  by 

ao=l,  a1=  -1/4  . a2=  7/32,  a3=  -21/64 

a4=  0.  7143531  708  , a§  = - 2.  07l  282248  . a6  = 7.  557254769 

1 a?=  -33.  32008068.  ag=  16 1 . 5948751  . a<?  = -1019.  868939 

a1Q  = 6845.  104932 

7 

Then  by  Laplace  inversion  of  (22) 


- $ H>s,  j b.  v(3/2»L  0(v33/2)  ( 


where  bQ,  bj.  ...  blQ  are  given  by 

v _ * v _ nr  ^iTT  /4 


bQ  = 1 . bj  = JT 


/ 4 , b2  = - 7i/60 


b3  = -7fTel"/4/5l2,  b4*-0. 


4398134  x 10‘ 


b5  = -0.  4109687  x 10‘3/T  e1*/4  . bg=  0.  1 1 2286 1 x 10' 3 i 


-5  r-  Jtt/4 


b?=  0.9182121  xlO  -fiT  e 


, bg=  0.  2093046  x 10 


(23a) 


(24a) 


b9  = 0.  I637812xl0"6/[r  e'llT/4  , bjQ=  -0.  3623427 x 1 0*7  i 
All  the  terms  9hown  have  been  used  for  subsequent  numerical  evaluation  of  G. 
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Actually,  the  inversion  yields  for  the  factor  outside  the  braces  in  (24)  the 
large  argument  approximation  of  (i/2)H^  *(ks) . the  Green's  function  for  an  infinite 
perfectly  conducting  plane.  By  inserting  the  exact  limiting  value  of  a — oo  , 
the  formula  in  (24)  may  be  applied  as  well  in  the  near  field  of  the  source. 

To  trace  the  transition  from  concave  to  convex  curvature  when  the 
distance  s along  the  surface  remains  fixed,  we  allow  a to  change  contin- 
uously from  positive  to  negative  values  via  the  complex  excursion 
0<arga<ir.  To  keep  s = a 1 4»  -<t» 1 1 positive,  it. is  necessary  simultaneously 
to  have  arg  |4>  -4>  ' | =-  arg  a.  Moreover,  one  must  continuously  deform  the 
integration  path  in  the  complex  t-plane  to  keep  the  integral  convergent 
(exp(iyt)  oscillatory)  when  y is  allowed  to  become  complex  according  to  the 
rules  stated  above.  These  considerations  lead  to  a straight  line  path  along 
which  arg  t = - ir/3,  2 it/3,  when  arg  a = it  . Changing  variables  t = p exp(i2Tr/3). 

one  finds 


G ~ 


iks  , 2/3  w (p) 

e (!«.  f _!_ 

2irka  2 w|(p) 

U 


^du 


(25) 


where  y is  again  given  by  (16a)  and  the  path  C proceeds  along  the  real  p 

axis.  The  expression  in  (25)  is  the  known  result  for  the  field  on  the  surface 

of  a convex  perfectly  conducting  cylinder.  By  contour  deformation  about  the 

singularities  p at  w',  (p  ) = 0,  one  derives  the  creeping  wave  series,  and  by 

pip 

expansion  analogous  to  (23)  and  Laplace  inversion  the  limiting  transition  as 
y — 0.  Thus,  one  may  track  the  field  continuously  as  the  curvature  changes 
from  concave  to  convex  between  fixed  source  and  observation  points  on  the 


surface. 
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III.  Numerical  Results 

Extensive  numerical  computations  have  been  performed  for  ka  = 10. 

30,  50,  100  and  1000  to  chock  the  accuracy  and  range  of  validity  of  the 
various  formulations  in  Section  II.  Results  are  presented  in  Figs.  6-12. 
for  ka  = 10,  100  and  1000  only  sin>.e  all  of  the  relevant  formulas  apply  to 
ka  = 30  and  50  with  accuracy  similar  to  th3t  obtained  for  the  larger  values 
ka  = 100,  1000. 

Except  for  the  low  (ka)  values  (ka  = 10),  the  whispering  gallery  mode 
and  continuous  spectrum  representation  in  (3),  for  arbitrary  y,  has  been 
taken  as  the  reference  solution,  with  which  the  other  results  are  compared. 
The  failure  of  this  representation  for  low  (ka).  as  evidenced  by  the  oscillations 
in  the  curves  in  Fig.  8(a)  for  small  y.  Is  attributed  to  difficulties  in 
evaluating  the  eigenvalues  of  the  whispering  gallery  modes  with  sufficiently 
high  precision;  for  the  larger  (ka)  values,  the  uniform  approximations  in 
(6)  are  adequate.  A remarkable  feature  of  this  representation  is  the 
sensitivity  of  the  field  to  the  exact  number  of  whispering  gallery  modes, 

as  evidenced  in  Fig.  6 for  ka  = 100.  Here.  = 32,  and  any  deviation  from 

this  number  causes  oscillations  about  the  true  field  magnitude. 

It  can  also  be  seen  from  Fig.  6 that  as  the  number  of  whispering 
gallery  modes  approaches  = 32,  the  field  magnitude  also  approaches  the 
true  value.  The  continuous  spectrum  is  necessary  to  establish  the  near 
field  for  moderate  y.  Analytically,  this  means  that  the  Neumann  function 
in  (4b)  together  with  the  discrete  spectrum  (wh?  oering  gallery  modes) 
produces  the  proper  singularity  near  the  source  point.  Omission  of  the 
continuous  spectrum  leads  to  errors  that  can  be  assessed  in  Fig.  7. 
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The  ray  plus  canonical  integral  representation  in  (14)  is  seen  to  pro- 
vide an  excellent  approximation  for  all  cases  (i.  e.  . for  ka  = 10,  100,  1000) 
as  shown  in  Figs.  8(a)-(c)  provided  that  the  number  of  rays  (N  + l)  is  chosen 
according  to  the  criterion  in  (16b).  Because  of  the  overlap  of  the  curves  for 
various  N.  no  difficulty  arises  in  switching  from  one  formulation  to  the  other 
as  y varies.  Tabulations  of  the  canonical  integral  in  (18a)  are  given  in 
reference  9 for  the  perfectly  conducting  boundary  and  for  non-vanishing  sur- 
face impedance.  When  ka  is  small,  this  representation  may  serve  as  a 
reference  because  it  does  not  exhibit  the  sensitivity  (see  Fig.  8(a))  encountered 
with  the  whispering  gallery  mode  and  continuous  spectrum  formulation.  As 
shown  in  Fig.  11(a).  the  ray  plus  canonical  integral  representation  for  small 
ka  agrees  well  with  the  near  field  expansion  in  (24). 

The  most  intriguing  and  physically  most  appealing  formulation  is  based 
on  the  mix  of  (N+l)  rays  and  M whispering  gallery  modes  in  (21).  The 
whispering  gallery  modes  plus  the  last  term  in  (21)  account  for  the  field  near 
the  surface  where  geometrical  optics  becomes  invalid.  Because  of  the  in- 
equality (21a)  that  must  be  satisfied  by  the  number  of  rays,  N+l.  the  inclusion 
of  multiply  reflected  rays  requires  larger  separation  v between  source  and 
observation  points.  For  any  given  N,  as  y increases  in  its  applicable  range, 
a new  whispering  gallery  mode  must  be  included  whenever  the  N-times  re- 
fleeted  ray  touches  the  caustic  of  that  mode.  The  process  for  N=0  is  depicted 
in  Fig.  9a.  The  concentric  circles  show  the  modal  caustics  of  the  whispering 
gall  ery  modes  while  the  angle  w^  represents  the  eigenvalue  of  the  K-th  mode. 
Within  the  ranges  of  P^,  P 2 and  P^.  one  must  include  the  M=1 . M=4  and 
M = K modes.  Corresponding  curves  for  |g  | are  shown  in  Fig.  10.  It 

should  be  noted  that  the  direct  ray  plus  whispering  gallery  mode  rep- 

I I 

resentation  yields  a good  approximation  for  y > 2(2JA)*  . 


This 
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representation,  however,  is  not  recommended  for  large  values  of  y since 
many  modes  must  then  be  included.  For  example,  N = 0,  M = 24,  at 
Y = 10. 

Inclusion  of  reflected  rays,  when  appropriate,  reduces  the  number  of 
required  whispering  gallery  modes,  as  depicted  in  Pig.  9(b).  Thus,  at 
the  observation  point  P(a,4>),  inclusion  of  only  the  direct  ray  (N=0)  requires 
six  whispering  gallery  modes  while  inclusion  of  the  direct  and  singly  re- 
flected rays  (N  = 0,1)  requires  only  a single  mode,  etc.  Fig. 9(b)  also  shows 
that  the  applicable  range  of  the  single  mode  plus  ray  representation  becomes 
wider  as  the  number  of  rays  increases.  As  seen  from  Pigs.  ll(b)-(c),  it 
is  possible  to  represent  the  field  accurately  with  only  the  lowest  (most 
tightly  bound) whispering  gallery  mode  provided  that  a higher  order  reflected 
ray  is  included,  when  necessary,  as  y increases.  Accordingly,  although 
(21)  is  a good  representation  for  the  ranges  satisfying  the  inequality  (21a), 
the'simplest  and  physically  most  appealing  expression  will  be  (21)  with  the 
inequality  ( 1 6b).  Figs.  ll(a)-(c)  have  been  obtained  with  the  inequality  (16b). 
These  results  are  surprisingly  simple  compared  with  the  whispering  gallery 
mode  and  continuous  spectrum  representation,  for  which  32  modes  for 
ka  = 100  and  319  modes  for  ka  = 1000  must  be  included  in  order  to  obtain 
accurate  field  values.  In  Figs.  ll(a)-(c),  the  number  indicated  along  the 
curves  should  be  read  as  follows:  for  example,  M=2  fN=l  ] .denotes  the 
range  wherein  the  direct  and  singly  reflected  rays  plus  2 whispering  gallery 
modes  are  applicable,  while  the  circles  indicate  the  starting  points  of  the 
relevant  intervals. 

Analytically,  the  emergence  of  a whispering  gallery  mode  (contact  of 
ray  and  modal  caustic)  is  characterized  by  the  coalescence  of  a pole  and  a 
saddle  point  in  the  remainder  integral  Therefore,  the  field  in  (21) 
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ia  discontinuous  at  those  points.  The  discontinuity  can  be  avoided  on  use 
of  a Fresnel  integral  transition  function  in  the  vicinity  of  these  \ values. 
However,  we  have  found  that  continuity  can  be  established  by  using  the 
simple  result  when  pole  and  saddle  point  coincide  (this  corresponds 
essentially  to  the  midpoint  between  the  two  discontinuous  endpoints)  and 
then  resorting  to  a perturbation  expansion. 

The  ray  and  whispering  gallery  mode  representation  is  invalid  in  the 
near  field  but  it  can  there  be  joined  to  the  near  field  expansion  in  (24).  Con- 
cerning the  latter,  it  should  be  noted  in  Fig.  12  that  inclusion  of  more 

perturbation  terms  docs  not  lead  to  continuing  improvement  or  extension  of 

21/2 

the  range  of  validity;  for  ka  = 100,  the  perturbation  terms  up  to  y seem 
to  give  the  best  result.  The  near  field  forms  with  8 perturbation  terms  are 
shown  in  Figs.  1 1 (a) - (c)  for  ka  = 10,  100  and  1000,  respectively.  From 
each  figure,  one  can  see  that  the  changeover  from  the  near  field  form  to 
the  direct  ray  plus  lowest  whispering  gallery  mode  (i.  e.  , M=l)  representa- 
tion can  be  performed  at  the  crossover  with  negligible  error. 
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IV.  Conclusions 

This  analytical  study  of  the  characteristics  of  the  surface  field  induced 
on  the  concave  side  of  a perfectly  conducting  circular  cylinder  by  a magnetic 
line  source  located  on  the  surface  has  furnished  a basic  understanding  oi  the 
propagation  mechanisms  that  are  operative  in  various  parameter  ranges 
involving  source  and  observation  point  location,  cylinder  radius  and  wave 
frequency.  Alternative  representations  of  the  field  solution  have  been  given 
and  their  relative  accuracies  assessed  by  extensive  numerical  calculation. 
'The  results  have  been  discussed  in  detail  in  Section  III.  They  permit,  in 
particular,  a quantitative  appraisal  of  when,  and  under  what  circumstances, 
geometrical  optics  begins  to  play  a role,  and  how  the  geometric  optical  field 
must  be  modified  to  account  either  for  near  field  effects  or  for  improperly 
formulated  ray  fields  that  have  undergone  many  reflections.  Quite  re- 
markably, the  latter  can  be  accounted  for  by  inclusion  of  only  a single  (the 
most  tightly  bound)  whispering  gallery  mode  and  a partial  ray  field,  although 
other  possibilities  exist  that  incorporate  several  whispering  gallery  modes. 
The  asymptotic  field  solutions  have  been  cast  in  a form  that  permits 

4 

generalization  to  variable  surface  curvature  or  variable  surface  impedance  , 
although  some  questions  remain  to  be  resolved. 

Calculations  for  surfaces  with  finite  surface  impedence  Z have  also 
been  performed,  as  have  tabulations  of  the  numerically  evaluated  canonical 

integral,  for  Z = 0 and  selected  values  Z 4 0.  These  results  will  be  _ 

9 

presented  elsewhere  . 


-20- 


Appendix  - The  Remainder  Integral 


As  noted  in  the  text,  the  remainder  integral  R^^  is  the  same  as 
in  (15)  except  for  the  replacement  of  C by  in  Fig.  4.  With  the 
aid  of  the  Wronskian  relation 


1 t /Jv<k*)  Hv(2)(ka)  \ 
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(26) 


Rj^N  can  be  represented  by 
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v ' \ N+l  iv 

»(2,(ka)/r  6 


dv 


(27) 


Along  the  contour  , one  may  use  in  the  integrand  the  Debye  asymptotic 

formulas  (see  (8),  and  similarly  for  the  undifferentiated  cylinder  functions) 
to  obtain 

RMN~  ^Ztt ^ tan  | ka[cos  w - (-|--w)  sin  w]+  -i^  elkaqN*w*dw 


(28) 


where  q^(w)  and  wg^  are  given  in  (10b)  and  (11),  respectively.  If 

(w^+^  + w]y£^/2,  ibe  asymptotic  evaluation  involving  the  tangent 
function  in  the  integrand  of  (28)  yields  a negligible  contribution  compared 
with  the  remaining  term.  In  this  case,  referring  to  (10a)  and  (12), 


• (29) 
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FIGURE  CAPTION’S 
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FIG.  1 
FIG.  2 
FIG.  2(a) 

FIG.  2(b) 
FIG.  3 
FIG.  4 


FIG.  5 


FIG.  6 
FIG.  6(a) 
FIG.  6(b) 
FIG.  7 
FIG.  8 


FIG.  8(a) 
FIG.  8(b) 
FIG.  8(c) 
FIG.  9 


Direct  and  multiply  reflected  rays  for  a circular  boundary. 
Reflected  ray  caustics. 

Caustics  A for  rays  reflected  n times;  these  ray  con- 
tributions are  confined  to  the  region  between  A and  the 
boundary. 

Caustic  formation  due  to  singly  and  doubly  reflected  rays. 
Whispering  gallery  mode:  modal  field,  modal  rays  and  caustic 


Integration  path  and  singularities  in  complex  v'  -plane. 

x --  zeros  v of  H1  ' 7 (ka ) : • • --  zeros  v of  J'  (ka). 

p v m v 

Cq  is  the  steepest  descent  path  through  v’sn* 

Integration  paths  and  singularities  in  complex  vv-plane 
(v  - ka  sin  w).  C1  corresponds  to  C in  Fig.  4 and  can  bo 
deformed  into  the  steepest  descent  path  C^  through  wg  . 

..  --  zeros  w of  J1  (ka);  the  polos  w are  not  shown, 
m v ' r p 

Whispering  gallery  (W.G.)  modes  plus  continuous  spectrum. 


Few  modes. 


Many  modes. 


Effect  of  continuous  spectrum 

Ray  plus  canonical  integral  representation.  Both  tin*  y and 
l0-0‘  I coordinates  are  indicated.  The  ranges  in  y and 
| | corresponding  to  N - - t (no  geometric -optical  ray) 

in  (14),  N = 0(1  ray),  N = 1 (2  rays)  and  N = 2(3  rays)  are  as 
shown. 


whispering  gallery  modes  and  continuous  spectrum 


JG 1 _ |R  1 

- I r\  4.  D 

1 U | - | j | 

* 

,gn=o+  rn=o  1 

ka  = 10 
ka  = 100 
ka  = 1000 

Geometrical  interpretation  of  direct  ray  plus  whispering 
gallery  mode  representation 


? 


FIG.  9(a)  N = 0 (direct  ray  only) 
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FIG.  9(b)  N = 0.  1,  2 

FIG.  10  Direct  ray  plus  whispering  gallery  mode  representation. 

M = K indicates  point  of  emergence  of  K-th  W.  G.  mode. 
The  solid  curve  is  calculated  from  (3). 


FIG.  11  Ray  plus  whispering  gallery  mode  representation.  The  solid 

curves  serve  as  reference  solutions.  Here,  the  |<f>-cf>'| 
coordinate  is  not  shown  since  it  has  been  depicted  in  Figure  8 


[Ns0]  IdJ'c  r c I 

PI-  |sGn,0  Gm| 


M =K  [ N = 1 ] 
o~ |g|  = 

M=1[N=2] 

|gI  s 


1 K 

n=0VJnTm=l  ~m"  2~N=1 


Z«G_  + S_JG  -iG, 


n“=0Gn+  Gm"  2°N=2 


near  field  form 


FIG.  11(a) 

FIG.  11(b) 
FIG.  11(c) 
FIG.  12 


ka  = 10  (M^=3).  The  solid  curve  is  calculated  from  (14)  and 

(18). 

ka  = 100  (M^=32).  The  solid  curve  is  calculated  from  (3). 
ka  = 1000  (Mj=319).  The  solid  curve  is  calculated  from  (3). 
Near  field  form,  ka  = 100 


FIG.  6 (a) 


FIG. 
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2.  High  Frequency  Fields  Excited  by  a Line  Source 
Located  on  a Concave  Cylindrical  Impedance  Surface 

by 

T.  Ishihara  and  L.  B.  Felsen 
Department  of  Electrical  Engineering 
Polytechnic  Institute  of  New  York 
Brooklyn.  New  York  11201 

Abstract 

A previous  study  of  high  frequency  currents  induced  by  a line  source 
on  a perfectly  conducting  concave  cylindrical  surface  is  extended  to  the  case 

of  non-vanishing  surface  impedance  Z . Alternative  field  representations 

s 

are  formulated  and  evaluated  asymptotically  as  combinations  of  ray  optical, 
whispering  gallery  mode,  surface' wave,  continuous  spectrum,  and  canonical 
integral  contributions.  Numerical  calculations  provide  an  insight  into  the 
accuracy  and  utility  of  the  various  formulations.  Sufficiently  far  from  the 
source  point,  a combination  of  ray  optical  fields  and  tightly  bound  whisper- 
ing gallery  modes  was  previously  found  to  be  a most  appealing  form  when 
Z =0.  As  the  surface  impedance  becomes  more  dissipative,  the  whisper- 

3 

ing  gallery  modes  are  weakened  by  attenuation  and  eventually  render  the 
ray  optical  fields  adequate  by  themselves.  A representation  in  terms  of 
rays  and  a canonical  integral  is  found  to  be  useful  for  all  parameter  ranges. 
The  canonical  integral  has  been  evaluated  numerically  and  tabulated. 


This  work  was  supported  by  the  U.S.  Army  Research  Office.  Durham. 
North  Carolina,  under  Contract  No.  DAHC  04-75 -G -0152. 


I.  Introduction 


When  high-frequency  fields  impinge  on  a concave  surface  with  large 
radius  of  curvature,  the  induced  fields  thereon  differ  markedly  from  those 
for  the  well  explored  convex  case.  A recent  study*  has  dealt  in  detail  with 
the  two-dimensional  problem  posed  by  line  source  excitation  on  the  interior 
boundary  of  a perfectly  conducting  circular  cylinder.  Alternative  field 
representations,  accounting  only  for  propagation  phenomena  on  the  surface 
segment  between  source  and  observation  points,  were  formulated  and 
evaluated  asymptotically.  The  field  expressions  so  obtained  were  inter- 
pretable physically  in  terms  of  constituents  involving  whispering  gallery 
modes,  ray-optical  fields,  near  field  effects,  continuous  spectrum  and 
canonical  integral  contributions.  Extensive  numerical  comparisons  showed 
the  accuracy  and  range  of  applicability  of  each.  This  investigation  has 
provided  fundamental  physical  and  quantitative  insight  into  the  field  be- 
havior on  the  perfectly  conducting  boundary. 

In  the  present  paper,  the  analysis  is  extended  to  accommodate  sur- 
faces with  non-vanishing  surface  impedance  Z . Because  the  alternative 

field  representations  (Section  II)  for  Z ^ 0 are  derived  by  the  same 

s 

techniques  as  those  for  Zg  = 0,  the  presentation  is  kept  concise,  with 
frequent  reference  made  to  the  earlier  study*  for  some  of  the  details. 
Representative  numerical  calculations  in  Section  III  show  how  increasing 
surface  impedance  affects  the  surface  field.  When  losses  are  appreciable, 
the  closely  bound  whispering  gallery  modes  are  rapidly  attenuated,  leaving 
a properly  formulated  ray-optical  field  as  the  dominant  and  adequate 


contributor. 
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II.  Alternative  Representations 

A.  Basic  Green's  Function 

The  relevant  magnetic  line  source  Green's  function  fc!r  the  interior 
of  a circular  cylinder  with  radius  a_  and  constant  surface  impedance  Z 
is  formulated  in  an  infinite  angular  space  wherein  the  azimuthal  coordinate 
4>  ranges  from  -ccto  +00.  2 With  the  source  point  located  at  (p*»  $' ) 
in  a cylindrical  (p.  <J>)  coordinate  system,  the  Green's  function  satisfies  the 
inhomogeneous  wave  equation 


(-,-pf  + ~ ^7  + k2)  G(p,  p ')  = - -6«,-<.’)6(p-p')  (i) 

p aP  p2  a<p2  - - p 

with  an  "angular  radiation  conditiod'  in  the  4*  direction,  and  the  radial 
boundary  condition 


= ikZ'G  at  p = 


a 


G finite  at  p = 0 (la) 


Here  Z'  = Z /Z  is  a normalized  surface  impedcnce,  and  Zq  is  the 

6 O 

impedance  of  free  space.  A time  factor  cxp(-iL't)  is  suppressed. 

Expressing  the  Green's  function  G as  a contour  integral  and  re- 
moving diffraction  effects  at  the  origin  arising  from  the  "angularly  matched" 
condition  , one  obtains  a modified  form  that  contains  only  the  essential 
propagation  phenomena  from  the  source  point  to  the  observation  point  in 
the  presence  of  a cylindrical  boundary  segment.  When  the  source  and 
observation  points  are  both  located  on  the  impedance  boundary  (i.  e.  , 
p = p'  = a),  the  modified  Green's  function  G becomes 

G=  1 f - exp(iv_U-»'[)  

i(Trka)2  C [j;(ka)-iZ'Jv(ka)][IF<2,(ka)-iZ'Hv(2)(ka)J 


dv  (2) 
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where  k is  the  free-space  wavenumber,  and  the  prime  on  the  cylinder 
functions  denotes  the  derivative  with  respect  to  the  argument.  The  contour 
C and  the  singularities  of  the  integrand  in  the  complex  v -plane  are  shown 
in  Fig.  1.  In  the  lossless  limit  Z'  * 0,  this  formulation  reduces  to  the  one 
in  reference  1« 

Anticipating  subsequent  contour  deformation  in  the  complex  v-plane  , 
we  examine  the  location  of  the  pole  singularities  determined  by  the 

resonance  equation 

J'  (ka)  -iZ'Jy  (ka)  = 0 , m = 1 , 2 M(1  Revm>  0 

vm  m 

These  poles  in  the  integrand  of  (2)  describe  whispering  gallery  modes.  For 

the  special  cases  arg  Z '=  + 90°,  the  roots  are  real.  If  for  arg  Z '=-90°, 

the  magnitude  |Z'  J of  the  normalized  impedance  is  greater  than  a certain 

value,  one  finds  for  the  first  root  m=l  that  v^>  ka.  The  corresponding 

field  contribution  then  describes  a surface  wav,e  that  can  exist  also  on  a 

plane  boundary.  It  is  interesting  to  note  that  all  of  the  whispering  gallery 

modes  represent  fast  waves  when  the  boundary  is  perfectly  conducting.  * 

However,  for  Z'  =£  0,  there  may  exist  slow-wave  type  whispering  gallery 

modes  with  Re  v > ka,  which  become  surface  waves  when  Im  v =0,  in 
m m 

addition. 


To  facilitate  the  numerical  evaluation  of  the  poles  in  (3),  it  is 
convenient  to  separate  the  domain  Re  >0  into  two  parts.  In  one  of 
these,  where  |y  |=  ka,  the  Fock  approximations  for  the  bessel 


functions  are  applicable: 


Jv(ka)~  (2/ka)  3Ai(t)  . 


H<1,Z\ka)  ~ (2/ka)  ~ w^  2(t)  . H^1,fi,(ka)  ~ -(2/ka)  ’ ^ 2(t) 


• Jj(ka)  ~ -(2/ka)  JAl' (t) 

ui  (1.2) 

IV  r,  \ ~ / "J  /V^\  3'  , 


where 


Wj  2(t)  = Ai(t)  + i Bi(t) , 


v = ka  + (ka/2)  *t 


r 


The  resonance  equation  for  the  first  M whispering  gallery  modes  bound 
most  closely  to  the  impedance  surface  is  then  approximated  by 

, i 

Ai'(t  ) + iaAi(t ) = 0,  m = 1 , 2,  • . • , M,  a = Z ( -ft  ' (5a) 


m 


m 


where  the  t are  solutions  of  the  differential  equation 
m 


dt 


m 


1 


" t -e2 

m 


| = -it* 


(5b) 


-=• 

In  the  other  domain,  where  |v  -ka| >0(  | v | 3 ) , one  may  employ  the  Debye 
approximations 


Jv (ka) 


- V—1 

A/  it  ka  c 


os  w 


sin  [4  (w)  + tt  /4] 


J^(ka) 


(1) 


- rjs 

■yj  ir 


cos  w 


ka 


cos  [£  (w)  + it  /4] 


(6) 


H ^(ka)  ~ /— t exp[+iC,(w)  + iir/4] 

v ' ' -v  irkacos  w 


(1) 

H,(Z)  (ka) 

v 


-ft 


cos  w 


ka 


exp  [+  i C(w)  + i^/4] 


•with 


;(w)  =ka[cosw-(ir/2-  w)sin  w],  v = ka  sin  w.  Re  w > 0 (6a) 

to  obtain  the  resonance  equation 


cos  w 


m 


cos[t;(w  ) + it  / 4 ] - iZ1  sin[£,(wm)  + tt  /4]  = 0,  m=M+l . . . . V) 


It  is  found  that  the  two  approximations  in  (5a)  and  (7)  have  an  operlapping 
domain  wherein  one  may  switch  from  one  to  the  other.  For  the  case  in 
Fig.  2,  M = 4 has  been  found  satisfactory. 


B.  Whispering  Gallery  (W.  G.  ) Mode  and  Continuous  Spectrum 
Representation 

When  the  contour  C in  Fig.  1 is  deformed  into  a contour  ex- 
tending along  the  imaginary  v-axis.  the  Green's  function  can  be  represented 
as  a sum  of  whispering  gallery  (W.  G.)  modes  Gm  and  a continuous  spectrum 


Mi 

2 . Gm  + RM. 

m=l  1 


(8) 


where 


Jv  (ka)  exp[ivm|4>-4>' l + iir/2]’ 
m 

kaO/av')7r  (ka)-  iZ'Jy  (ka)} 


with  the  satisfying  the  resonance  equation  in  (3).  In  obtaining  the 

residue  contributions  in  (9)  due  to  the  poles  at  , the  Wronskian 
relation 


1_ 1 r Jv  fra)  H(2)(ka) 

iirka  [ H^(  3 (ka)  -iZ ' H j 3 (ka)  ] [ r (ka)  - iZ ' J (ka)  ] 2 [ F (ka)  - iZ 1 J (ka)  H^(ka)  - iZ ' H®(ka) 

(9a) 

has  been  utilized*  Using  the  Fock  approximations  in  (4)  and  (5),  one 
finds 


exp(iks+ivt  -itr/2)  

G ~ r — — j . m ■ i . • • M 

m 2(ka/2)  J(tm  + a) 


(10) 
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r 


whe  re 


,ka  s 
7 = (“T  ) . T 


(10a) 


Here,  s denotes  the  path  length  between  source  and  observation  point 
along  the  boundary,  and  y is  an  arc  length  parameter.  For  the  remaining 
modes,  (6)  and  (7)  may  be  used  to  obtain 


exp(ikssinw  +itr/2) 

- - - ’ ' ~~  **2 

ka(ir  /2-w  )cos  wm[  1 -(Z'/cos  wm)  ) 


m-M  + 1 . . . M, 


From  the  behavior  of  the  roots  as  discussed  in  Section  A,  one  can  see  that  the 
residue  contribution  arising  from  the  m=l  pole  singularity  at  tj  (the  first 
root  of  (5))  describes  either  a whispering  gallery  mode  bound  most  closely 
to  the  impedance  boundary  or  a surface  wave  mode  when  Re  tj  > 0, 

Im  tj  = 0.  In  the  latter  instance,  in  (10)  becomes  a surface  wave  field. 

As  in  the  lossless  case,  the  continuous  spectrum 


Rm  = -H- 7 -dv  (12) 

M1  i(trkaj  -iso  [JJ  (kaJ-iZ'J^  (ka)]  [H^  (ka)-iZ'Hj'  '(ka)] 


can  be  approximated  by 


I w 

1 I 

ir  J 


exp( -ks  sinh  v) 


1 n o 1 - (Z.,/cosH\r)< 


(12a) 


C.  Ray -Optical  plus  Slow  Wave  Representation 

With  a view  toward  a ray-optical  representation,  it  is  convenient 
to  deform  the  contour  C in  Fig.  1 into  C,  whence  G becomes 


G = G + Gj  U(tt) 


A 


where  G is  the  same  as  G except  for  the  replacement  of  C by  C in  (2), 
and  Gj  is  the  contribution  (if  any)  from  the  pole  singularity  at  tj  located 
between  C and  C.  The  unit  step  function  U^)  contributes  only  when 
the  tj  pole  lies  between  the  two  contours.  The  extraction  of  Gj,  which 
may  be  either  a slow  wave  W.  G.  mode  or  a surface  wave,  is  desirable 
since  slow  wave  fields  do  not  exhibit  ray-optical  properties. 

To  convert  G into  ray-  optical  form,  we  substitute  the  traveling 

wave  expansion 


1 


JJ  <ka)-iZ'Jv  (ka) 


H^(2)(ka)-iZ'Hv(2)(ka)  n=0 


X (*r)n.  | r | < i (14) 


r = 


Hj(1)(ka)-iZ'H;(1)(ka) 

H'v(2)(ka)-iZ'Hv(2)(ka) 


(14a) 


into  (2),  with  G and  C replaced  by  G and  C,  respectively.  Then  with 

the  Debye  asymptotic  formulas  in  (6),  changing  variables  from  v to  w via 
v = ka  slnw  , 


G = 


n=0 


where 


c .fcH 

n ir 


n 2 f 7 

I r COS  W . I COS  W-Z 

1,  , , 7I>2  l cos  w+Z 

C'  (cos  W + Z')  \ 


-zA  ikacin(w) 


and 


q (w)  = | <b- 4>*  | sin  w + 2(n+l)[cos  w-(n/2-w)sinw] 


(15) 


dw  (15a) 


(15b) 


The  integral  in  (15a)  can  be  evaluated  in  terms  of  the  contributions  from 

the  saddle  points  w gn  and  wg  of  qn(w),  determined  from  the  condition  dqn/dw=0: 


w 


sn 


2(n+l ) 


— TT 

W = rr 
S 2 


(16) 
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A typical  path  passing  through  the  n-th  saddle  point  is  shown  in 
Fig.  3 (note  that  the  whispering  gallery  mode  poles  w have  been  eliminated 
from  the  integrand  in  (15  a)).  Upon  using  the  conventional  saddle  point 
formula  (ref.  2,  p.  382)  one  finds 


**P<‘kDnl 

(cos  w6n  + 


where 

D = 2(n+l ) a sin  [ | 1 /2(n+l )] 

n 


(17) 


(17a) 


and  T is  the  boundary  reflection  coefficient  for  ray  species  n. 
n 

cosw  - Z' 
r _ sn 

n “ cos  w +Z' 
sn 


(17b) 


The  expression  in  (17)  corresponds  to  a geometric -optical  ray  field  that 
has  undergone  n reflections  at  the  impedance  boundary  (see  Fig.  4).  The 
sum  in  (15)  can  actually  not  be  extended  to  n — oc  since  the  Debye  approxi- 
mations in  (6)  become  invalid  as  w -*  it  / 2.  This  situation  arises  either 
when  the  observation  point  approaches  the  source  point  so  that  J <}> 0 , 
or  when  the  number  of  rays  increases,  with  | <t> - ( fixed.  The  latter  case 
can  be  avoided  by  truncating  the  number  of  rays  at  some  n = N such  that 
wgN  is  sufficiently  less  than  tt/2  (see  Sec.  D)  while  the  former  case  will 
be  discussed  in  Section  F. 


D.  Ray  plus  Canonical  Integral  plus  Slow  Wave  Representation 
When  instead  of  (14)  the  partial  expansion 


Jj  (ka)  - iZ'Jy  (ka) 


H'v(2)(ka)-iZ'H  (2)(ka) 


f (-r)n  + jr 

n=0  v 


Izll 


N + l 


y (ka)  - iZ'J,  (ka) 


(18) 


is  substituted  into  (2)  (with  G and  C replaced  by  G and  C,  respectively), 
one  has  instead  of  (13)  and  (15) 


+ R 


N 


+ Gj  U(tj) 


(19) 


whe  re 

R _ 1 r (-r)N+1  exp(i  v | 6-<h'  |j dv 

N i(irka)2  l [H'v(2)(ka) -iZ'Hv(2)(ka)][J^(ka) -iZ'Jv(ka)] 

(19a) 

G , given  in  (17),  again  describes  a geometrical  ray  field  having  undergone 
n reflections.  The  remainder  term  accounts  for  ray  fields  reflected 

more  than  N times. 

One  way  of  dealing  with  (19a)  is  to  investigate  how  fast  the  integrand 
decays  in  the  neighborhood  of  the  point  v = kzt;  rapid  decay  permits  effective 
truncation  of  the  integration  path  and  the  use  of  further  approximations. 
Proceding  as  in  the  Z'  = 0 case,  * one  may  show  that  the  main  contribution 
to  the  integral  arises  from  the  neighborhood  of  v = ka  provided  that  the 
arc  length  parameter  y defined  in  (10a)  and  the  number  of  rays  (N+l) 

' satisfy  the  left-hand  portion  of  the  inequality 


■ - 1 <(N+1 ) < V-r  • (20) 

2(23A)*  2(2  3A)2 

where  A = A (ka,  Z')  = 0(1)  and  |z'|  «.  1.  The  right-hand  portion  of  the 
inequality  validates  use  of  the  Debye  approximations  in  the  first  term  on 
the  right-hand  side  of  (18),  which  has  been  used  to  derive  Gfi  in  (17).  The 
complete  inequality  in  (20)  provides  a criterion  for  the  number  of  geometric 
optical  ray  fields  that  are  identified  in  (19).  The  precise  choice  of  A is 
not  too  important;  a value  A » 2 has  been  found  adequate  in  our  various 
numerical  comparisons*  ( Sec.  III). 


Since  the  contributing  portion  of  the  integrand  in  (19a)  is  localized 


near  v - ka,  one  may  substitute  the  asymptotic  approximations  (-1)  to  obtain 

(21) 


2iir  ka 


Here,  l^j(\.a)  is  the  canonical  integral 


f ( w'j  (t)  + i«j\v((t)}/{\v^(t)  + iow,(t)}]N  + 1 


IN(v,a)  L [w'(t)  + iaw_(t)j  fAl'(t)  + inAi(t)J  e V dt  (22) 


with  the  path  in  the  complex  t plane  obtained  from  C in  the  complex 
v -plane  by  the  shift  implied  in  ( 4a  ).  The  canonical  integral  in  (22)  has 
been  evaluated  for  relevant  ranges  of  y and  three  values  of  a.  The  results 

5 

are  shown  in  Fig.  11  and  have  been  tabulated  . The  canonical  integral  for  a 
related  boundary  value  problem  corresponding  to  a different  integrand  and 
to  a = 0 has  been  tabulated  by  Babich  and  Buldyrev.  ^ 


E.  Ray  plus  Whispering  Gallery  (W.  G.  ) Mode  plus  Slow  Wave 
Respresentation 

An  alternative  way  of  treating  (19a)  is  to  deform  the  contour  C into 
in  Fig.  3 and  thereby  representing  R^as  a sum  of[  M-U(t  j )]  W.  G.  modes  whose 
poles  lie  between  C and  C^,  plus  some  other  remainder  term 
Referring  to  (9a),  one  finds 

RN^_Gm  + RMN  ’ ™=1+U(tl>  <23> 

m=m 


where 
RMN=  2 


.rka  / 


Jv  (ka) 


^(2)(ka) 


r (ka)-iZ'Jy  (ka)  H'(2,(ka)-iZ'H  (2)(ka) 
N [ v ' v 


(_r)N+leiv  U-4>' 


(24) 


The  sum  in  (23),  with  G given  in  (10)  or  (11),  describes  W.G.  mode 


dv 


-contributions 
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3 


The  evaluation  of  in  (24)  may  proceed  by  the  saddle  point 

technique,  after  replacement  of  the  cylinder  functions  by  their  Debye 
approximations  as  in  (6).  The  integral  then  becomes: 


..N+l  f 
•>)  f £° 

2 IT  1 CO 

CN  i 


oswtan{r.(w)  tir/4)  icosw  /cos  w - Z1 

osw-iZ1  tanjt,  (wj  + tt/4 } cos  w+Z'  l cos  w+Z' 


f \ N+l 

cosjv^Z'\  eikaqN(w)di 


where  f,  (w)  is  given  in  (6a)  and  q^(w)  in  (15b).  For  observation  points  such 
that  the  saddle  point  wsjyj  of  (25)  is  located  far  from  the  pole  singularities 
of  the  integrand,  and  for  a low  loss  impedance,  the  contribution  from  the 
first  term  inside  the  square  brackets  in  (25)  is  negligible  when  compared 
with  that  from  the  second  term.  Then  from  a comparison  of  (25)  with  (15a) 


and  (17) 


- 1 7? 


Thus 


■ if  S*1 


«.  *2 


cos  w 


P 1_  /—  /<  2_ 

m ' 2 N u " cos  w. 


m=m 


F.  Near  Field  Form 


As  observed  at  the  end  of  Section  C (see  also  (20)),  any  Green's 

function  representation  that  includes  geometric  optical  terms  fails  in  the  near 

I I 

field  where  the  arc  length  parameter  -y  in  (10a)  is  less  than  2(23  a)2.  We 
therefore  go  back  to  the  basic  integral  representation  in  (2).  Since  the  main 
contribution  to  the  integral  arises  from  the  vicinity  of  v - ka,  the  cylinder 
functions  may  be  replaced  by  their  Fock  approximations  in  (4)  . 

Then  using  the  Wronskian  relation  for  the  Airy  functions  (analogous  to  (9a))  and 

I 

applying  Cauchy's  theorem  in  the  upper  half  of  the  complex  t-plane,  one  finds 


-exp(iks)  ,ka  . * +*f  ^ Ai(t)exp(i\t) 

C 2rr  ka  2 _J_i6  Ai'(t)  +iaAi(t) 


Employing  large  argument  expansions  for  the  Airy  functions,  one  may 
generate  an  asymptotic  expansion  of  the  integrand  in  (28): 


Ai(t)  ~ ^ a± 

Ai'(t)  + aAi(t)  Aj  t(j+t)/2 


+ o(t'23^2) 


where  all  of  the  terms  retained  (i.  e.  , up  to  j = 21)  were  employed  in 
subsequent  numerical  calculations.  The  coefficients  aQ,  a^ a^ 


are  given  by 


IS  d (ia)3f  , j = 3i 

i=  0 J* 

k 3/  +1 

2 d.  . (iaj  * +1  , j = 3/+1 

1=0  -1'1 

A 3f+2 

E d.  .(ia)3Z+/:  . j = 31+2 

i =0  J* 


where 


dj.L=‘  • 

dj>L-l=-<j-2)/4  * 

d6.0  = 7/32 

dj . L-  2=dj  - 1 , L-  2 

+ 9/32  + 2 (j - 7) / 32,  for 

j = 3L+1.  3L+2,  (L> 

dj,L-2=dj-!.L-3 

+ 9/32  + 2(j -7)/ 32.  for 

j = 3L  (L  > 3) 

d,0= -21/64 

d10.0  = -49/64 

dll_0  = -85/64 

duo=0.  71435 

d12jl=  -130/64 

d,3  0 = 1. 64062 

d13l=-  185/64 

<3,4,0=  2-811982 

dul= -251/64 

dI5io  = -2-°7128 

d.  c .=  4.  29687 

13,  1 

di5i  2=  -329/64 

d,6  0= -4.  64329 

d16 .1  " 6- 1 2°59 

d,6.2=  -420/64 

dj7tos  _7’  81345 

d(7  8.  34373 

dl,_  2=  -525/64 

d18.0=  7-55725 

• d,8  ,=  11.  69283 

d,8  2=ll-°2293 

d18, 3=  "645/64 

d19.0=  ,6'  57034 

d,,  ,=  -1 6, 40t  7 

d19  2=  14.  21873 

dl9  3=  -.781/64 

d2 o.o- 27' 29129 

d2Q  1=  -22.09779  d20  2=  17.99556 

d20.3  = -93V64 

d21  o=_33,  32008  d21  1=40-0!068 

d21 , 2=  -28'  92270 

d2i  3=  22.  42181 

d21.4=  -1105/64 

7 

by  Laplace  inversion  of  (28),  with  (29). 

iks+irr/4  1 
G ~ e I 

i+l*  b </2  + olv11) 

1 . Ul<  1 

^ 2irk  s j 

where 


b.  = 

J 


J!t  a.  e^  4 ^ )!  . j odd 


►j/2a.  e1^^  / n (Em-1),  j 
3 / m=l 


(30) 


(30a) 


even 


The  restriction  |a|  < 1 in  (30)  arises  from  the  condition  | a Ai(t) /Ai’  (t)  | < 1 . 
required  for  validity  of  the  expansion  in  (29).  In  view  of  the  definition  of 
a in  (5a),  the  solution  (30)  can  therefore  not  be  applied  to  the  infinite 


plane  limit  a co  if  Z'  differs  from  zero. 
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III.  Numerical  Results 

Extensive  numerical  computations  have  been  performed  for  ka  = 100 
and  for  various  surface  impedance  values  to  show  how  the  amplitude  and 
phase  of  the  surface  impedance  affect  the  field,  and  to  check  the  accuracy 
and  range  of  validity  of  the  various  formulations  in  Section  II. 

Concerning  the  propagation  constants  for  the  W.G.  modes,  we  have 
obtained  solutions  of  the  differential  equation  (5b)  by  the  Runge-Kutta 
method.  The  results  are  shown  in  Fig.  5.  Here  tj,  t^  and  t^  are  the  first 
three  roots  of  Ai*(tm)  = 0 (i.  e.  , a = 0 in  (5a))  and  t^,  t 2 and  t 3 are  those 
of  Ai(tm>  = 0 (i.  e. , a = coin  (5a)).  The  solid  lines  and  the  dotted  lines  are 
equi-phase  and  equi-amplitude  curves,  respectively.  The  arrows  on  the 
equi-phase  curves  indicate  the  directions  along  the  root  loci  with  increasing 
amplitude  |||.  It  is  interesting  to  note  in  Fig.  5 that  the  first  zero  t 
crosses  the  imaginary  axis  for  some  values  of  £ , thereby  furnishing  the 
slow  wave  type  of  whispering  gallery  mode.  In  particular,  when  Re  t^>  0, 
Im  tj  = 0,  the  contribution  from  t^  describes  a surface  wave  mode. 

Except  for  arg  Z*  = + 90°  , the  imaginary  parts  of  the  whispering  gallery 
mode  poles  become  large  as  the  amplitude  of  Z'  increases  within  the 
ranges  under  consideration  (i.  e.  , jz'|  <0.  6).  It  should  be  noted  in 
Figs.  2 and  5 that,  for  given  Z'  and  arg  Z'  £ + 90°  , the  eigenvalues  for 
modes  bound  close  to  the  surface  have  large  imaginary  parts  compared  to 
those  located  far  from  the  surface,  thereby  de-emphasizing  the  importance 
of  the  most  closely  bound  whispering  gallery  modes.  This  is  to  be  expected 
on  physical  grounds  since  the  field  of  the  tightly  bound  modes  is  in  close 
proximity  to  the  lossy  guiding  surface.  It  can  also  be  seen  from  Fig.  2 
that  the  two  approximations  given  in  (5a)  and  (7)  have  an  overlapping  region 
wherein  one  can  switch  from  one  to  the  other.  M=  4 has  been  found 
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satisfactory  for  our  evaluation  of  the  field. 

In  Figs.  6-7,  we  have  compared  the  various  formulations  in 
Section  II  for  a purely  inductive  impedance  case  (|z'|  = 0.  22  and 
arg  Z'  = -90°)  to  check  accuracy  and  range  of  validity.  For  this  case,  the 
surface  wave  mode  ss  excited  since  the  first  root  of  (5a)  or  (5b)  is  located 
on  the  positive  real  axis  in  the  t-plane.  The  ray  plus  canonical  integral 
plus  slow  wave  representation  in  (19)  is  seen  to  provide  an  excellent 
approximation  provided  that  the  number  of  rays  (N  + l)  is  chosen  accoiding 
to  the  criterion  in  (20).  Because  of  the  overlap  of  the  curves  for  various 
N,  switching  from  one  formulation  to  the  other  can  be  performed  smoothly. 
The  ray  plus  canonical  integral  contribution  curves  are  shown  separately 
to  exhibit  the  effect  of  the  surface  wave  on  the  total  field. 

In  Fig.  7,  the  near  field  form  in  (30)  and  the  ray  plus  whispering 
gallery  mode  plus  slow  wave  representation  in  (27)  are  compared  with  the 
whispering  gallery  mode  (including  surface  wave)  plus  continuous  spectium 
representation  in  (8).  It  is  interesting  to  observe  that  for  large  enough  y , 
the  total  field  can  be  accurately  represented  just  by  the  surface  wave 
(M=l)  and  by  rays  provided  that  the  appropriate  number  of  reflected  rays 
is  included  as  y or  [ cp  — 4> * | increases. 

Curves  for  a high  loss  case  are  depicted  in  Figs.  8(a)  and  8(b), 
based  on  the  ray  plus  canonical  integral  representation  and  on  the  ray- 
optical  representation,  respectively,  with  the  latter  implying  omission 
of  the  canonical  integral  contribution.  Having  in  mind  an  application  to 
ground  wave  propagation,  the  high  loss  impedance  Z'  = 0.22  exp  (-i30.  5 ) 
has  been  obtained  by  assuming  a wet  ground  surface  with 
conductivity  o-  — 10  mho  m and  dielectric  constant  c ^ - 10  € ^ * 

subject  to  a wave  frequency  f = 10MHz.  Also  shown  in  Fig.  8(b)  is  the 
near  field  form  for  small  y . It  should  be  noted  that  in  (19)  docs  not 


contribute  since  no  pole  singularity  is  located  between  the  contours  C and 
C in  Fig.  1 (see  also  Section  II- A).  All  of  the  modes  here  are  fast  waves 
as  seen  from  Fig.  2.  Note  that  the  ray-optical  contribution  alone  is  adequate 
for  describing  the  field  sufficiently  far  from  the  source  point  because  of  the 
dissipative  attenuation  of  the  modes  bound  close  to  the  surface,  as  mentioned 
previously.  •' 

The  solid  curves  in  Figs.  6-8  are  obtained  from  (S)  with 
=32(ka=100).  This  form  of  the  solution  is  taken  as  a reference  for  all 
y . It  was  noted  in  the  study  of  the  perfect  conductor  case*  that  the 
numerical  values  derived  from  the  W.  G.  mode  plus  continuous  spectrum 
representation  are  extremely  sensitive  to  the  exact  number  of  modes  re- 
quired. Unless  all  32  modes  (for  ka=100)  are  included,  the  perfect  con- 
ductor curve  was  found  to  deviate  appreciably  from  the  correct  shape.  The 
same  remark  applies  to  the  finite  surface  impedance  case  in  Figs.  6-8 
since  the  higher  order  modes  have  small  attenuation  coefficients  (see  Fig.  2) 
even  though  the  losses  may  be  large.  Because  of  this  feature,  it  is  pre- 
ferable to  use  one  of  the  other  representations  for  field  calculation. 

In  Figs.  9(a)-(e),  we  have  shown,  for  given  phase  of  the  normalized 
impedance  Z',  how  increasing  the  amplitude  of  Z'  affects  the  surface 
field.  Curves  are  depicted  based  on  the  ray  plus  canonical  integral  plus 
slow  wave  (if  necessary)  representation,  and  on  the  near  field  form  for 
small  y . In  our  computations,  | Z * [ =•  0.  05  has  been  chosen  as  a sample  of 
a small  impedance,  and  |z'  J = 0.  22  as  a sample  of  a large  impedance.  If 
arg  Z*  = -90°  , corresponding  to  a purely  inductive  impedance  boundary,  the 
field  amplitude  increases  as  |z'|  increases  (see  Fig.  9(a))  while  the  con- 
verse is  true  when  the  boundary  impedance  is  purely  capacitive  (see 
Fig.  9 (e)).  This  behavior  may  be  attributed  to  the  greater  and  lesser 
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confinement,  respectively,  of  the  lowest  order  modal  field  in  the  two  cases. 
As  the  imaginary  part  of  roots  increases  (see  Fig.  4), the  field  transmitted 
to  an  observation  point  tends  to  decrease  (Figs.  9(b)-(d)).  Because  of  the 
absorption  on  the  surface,  the  interference  effects  between  rays  are  also 
weakened,  as  seen  from  the  decreasing  magnitude  of  oscillation  in  the 
curves  in  Figs.  9(b)-(d). 

/ 

The  curves  in  Fig.  10  show  the  effect  of  the  phase  of  Z'  on  the 
field  when  the  amplitude  remains  constant  at  |z'|  = 0.05.  One  observes 
that  the  field  tends  to  increase  as  the  phase  of  Z'  deviates  from  zero;  this 
could  have  been  predicted  from  the  behavior  of  poles  in  Fig.  5.  Note  that 
the  maxima  and  minima  of  the  field  shift  backward  for  negative  phase  (i.  e. , 
-90°  < arg  Z'  < 0°)  and  forward  for  positive  phase  (i.  e. , 0°  < arg  Z1  <_  90°). 


IV.  Conclusions 


A variety  of  alternative  representations  have  been  explored  in  this 
study  of  the  surface  fields  excited  by  a line  source  located  on  a concave 
cylindrical  impedance  boundary.  As  for  the  previously  investigated 
perfectly  conducting  case',  these  different  formulations  provide  new  in- 
sights into  the  propagation  mechanisms  that  prevail  for  various  parameter 
ranges  involving  the  boundary  shape,  separation  of  source  and  observation 
points,  and  the  surface  impedance.  From  the  detailed  discussion  and  in- 
terpretation of  the  numerical  results  in  Section  III  emerges  the  role  played 
by  the  geometric  optical  field  and,  for  sufficiently  inductive  boundaries,  by 
the  surface  wave.  The  presence  of  the  surface  wave  leads  to  an  enhance- 
ment of  the  field  over  that  observed  on  a perfect  conductor.  When  losses 
are  appreciable,  attenuation  of  mode  fields  bound  close  to  the  boundary 
establishes  the  geometric  optical  field  as  the  dominant  and  adequate  con- 
stituent. The  versatility  of  representations  involving  the  canonical  integral 
has  again  been  confirmed,  and  tabulations  have  been  provided  for  several 
values  of  surface  impedance. 
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1 Integration  path  and  singularities  in  complex  v -plane. 

xx'-zeros  v of  [H1^  \ka)  - iZ'Il^  ^(ka)];  »■»  --  zeros  v of 
p v V na 

I r (ka)  - iZ'Jy  (ka)]. 

2 Zeros  of  [J’v  (ka)-iZ'Jv  (ka)]  in  complex  w-^plane. 

• •--zeros  for  Z ' = 0 (perfect  conductor);  oo- - zeros  for  Z '=  0.05  exp(- 
(small  loss);  3 3 , AA  --  zeros  for  Z'  = 0.  22  exp(-i30.  5°) 

(large  loss).  The  first  five  roots  AA  --  are  obtained  from  (5b) 
and  compared  with  0 Q --  obtained  from  (7). 

3 Disposition  of  integration  paths 

3(a)  Steepest  descent  paths  in  complex  v -plane 

3(b)  Steepest  descent  paths  in  complex  w-plane  (v  =ka  sin  w) 

4 Direct  and  multiply  reflected  rays  for  a circular  boundary. 

Geometrical  quantities  Dn,  D , D. ,,  w , w and  w , are  also 

UlcsOsl  s<- 

shown. 

5 Roots  of  jAi'(t  ) + iaAi(tm)]  as  obtained  from  (5b).  The 
numbers  on  the  solid  and  dashed  curves  denote  the  phase  and 
amplitude  of  £ = -i  a , respectively. 

6 Ray  plus  canonical  integral  plus  slow  wave  representation  in 
(19),  with  Z'=  0:  22  exp(-i90°).  Both  the  \ and  |ct>-<J>'  | 
coordinates  are  indicated.  Also  shown  are  the  ranges  in  y and 

J 1 corresponding  to  N = -1  ( ■■■«■  c>-  ) (no  geometric -optical 

ray);  N = 0 ( ) (1  ray);  N = 1 ( — • — )(2  rays);  N = 2 ( — ) 

(3  rays).  The  solid  curve  is  calculated  from  (8)  and  serves  as 
the  reference  solution.  The  heavy  curves  include  the  surface 
wave  in  (19)  (U(tt)  = 1 for  this  case),  while  the  light  curves  have 
been  obtained  from  (19)  with  the  surface  wave  omitted. 
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FIG.  8(a) 
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FIG.  9 


Ray  plus  IV.  G.  mode  plus  slow  wave  representation  in  (27), 
and  near  field  form  in  (30)  (--),  with  Z'=  0.  22  exp(-i90°). 

The  reference  solution  (solid  curve)  is  calculated  from  (8). 

The  numbers  indicated  along  the  curves  should  be  read  as 
follows:  for  example,  M = 2 lN=l]  denotes  the  range  wherein 
the  direct  ray  and  singly  reflected  rays  plus  2 modes  (1  W.  G. 
mode  and  1 surface  wave  mode)  are  applicable,  while  the  circles 
indicate  the  starting  point  of  the  relevant  intervals.  The 

1 coordinate  is  not  shown  since  it  has  been  depicted  in 
Fig.  6. 

Ray  plus  canonical  integral  representation  in  (19),  and  near 
field  form  (--)  in  (30),  with  Z'  = 0.  22  cxp(-i30.  5°)  . Note  that 
U(tj)  = 0 in  this  example  (see  Fig.  2).  The  reference  solution 
(solid  curve)  is  calculated  from  (8). 


Canonical  integral  contribution  included, 

N=  0, N = 1;  N = 2. 


>-  N=-l ; 


Canonical  integral  contribution  R^  omitted.  — • • — N = 1; 
N = 2. 

Influence  of  the  amplitude  of  the  surface  impedance  Z',  cal» 
culated  from  (19)  and  (30)  (--),  for  various  phase  angles.  The 
heavy  curves  represent  a large  impedance  ( | Z 1 1 = 0. 22)  and  the  light 
curves  a small  impedance  ( | Z' | = 0.  05).  Here,  U(t  ^ ) = 0 except 
for  the  case  Z'=  0.  22  exp(-i90°).  — N = -1;  — — N = 0; 

— • • — N = 1;  and  — • • •—  N = 2. 
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Influence  of  the  phase  of  the  surface  impedance 

from  (19),  with  |z*|  =0.05. 
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FIG.  10(a)  Negative  phase  angle. 

arg  Z*  = 0°;  — — arg  Z'  = -45°; 

— • — arg  Z ' = -90°. 

FIG.  10(b)  Positive  phase  angle. 

arg  Z 1 = 0°;  — • — arg  Z'  = 45°; 

— • — •— argZ'  =90° 

FIG.  11  Plot  of  canonical  integral  ll^  (y  , or)  | vs.  y . 

Solid  curves:  . Z'  = 0 (a  = 0)  (perfect  conductor); 
- Z'  = 0.  05  exp(-i30.  5°)  or  Reo-  = 0.  16799 

lm  a = - 0.  09903  (small  loss)  ; r 

Zf  = 0.  22  exp(-i30.  5°)  or  Rea  = 0.  6997, 
lm  a = -0.  41261  (large  loss). 
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FIG.  9 (a) 


